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Abstract 

We study the class of neutrino mass matrices with a dominant block 
but unspecified 0(1) coefficients, and scan the possible models by the 
help of random number generators. We discuss which are the most 
common expectations in dependence of the adjustable parameter of 
the mass matrices, e, and emphasise an interesting sub-class of models 
that have large mixing angles for atmospheric and solar neutrinos, and 
an angle #13 close to the experimental limit. For those models where 
the lepton mass matrices are subject to Froggatt-Nielsen U(l) selection 
rules, we show that the neutrino mixing matrix receives important 
contributions from the rotations operating on charged lepton sector, 
which increase the predicted value of the angle #13 and the ee-entry of 
the neutrino mass. 



A specific, simple form of neutrino mass matrix has stimulated many studies 
in the last few years (see e.g. § |, 1 1, |, |, g | |, 0): 

M v <x diag[e, 1, 1] 0(1) diag[e, 1, 1] (1) 

where it is understood that e is a small parameter^], and 0(1) is the matrix 
of the "coefficients of order unity" (an argument for eq. ([j]) with a minimum 
of theoretical assumptions is presented in |l^|). The attempt to understand 
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1 The case e = 1 has been given theoretical support in jnj , and was discussed in detail 
in H and M; however, we regard it as an extreme case. 



the gross features of the fermion (neutrino) mass matrices by the help of 
factor e but without stressing too much the role of the coefficients is in the 
line of thought drawn by Froggatt and Nielsen p3| j. In their language, the 
matrix in eq. (|l|) is characterised by the presence of two degenerate U(l) 
charges, and this implies that a block of the mass matrix has relatively large 



entries ("dominant block"). Recently, the early suggestion [13] that the 
coefficients should be thought as random numbers was taken quite literally 
in the studies || [7|, ^, 1C]. We find this type of approach of interest, since 



it amounts to a scan of the possible theoretical models; the most frequent 
cases (in dependence on the value of e) are emphasised in this approach. We 
show in the appendix our assumptions on the coefficients, seen as random 
numbers, and discuss related problems. 

We begin by discussing the value of e putting emphasis on available neu- 
trino oscillation data. Postponing the interpretation of LSND indications, 
we know with good confidence that^ the atmospheric mixing angle #23 is in 
the range 45° ± 10°, while #13 belongs roughly to (0 — 10)°; the correspond- 
ing Amli is in the range (1.5 — 5) x 10 -3 eV 2 . (It is rather evident that the 
two informations on the mixing angles are in reasonable agreement with the 
mass matrix in eq. (|]), if e is sufficiently small). The situation with solar 
neutrino observations is less clear. We assume that the observations indicate 
again oscillations of massive neutrinos, with Am^ <C Am^, and select three 
regions of parameter space for further discussion: (1) a SMA region, namely 
the rectangle with vertices (4 x 1CT 6 ,2 x 1CT 4 ) and (1 x 10 _5 ,3 x 1CT 3 ) in 
the (Am^/eV 2 , tan 2 #12) plane, (2) a lma region, namely the rectangle with 
vertices (8 x 1CT 6 ,0.15) and (3 x 1CT 4 ,0.75) in the same plane, (3) a low 
(up to quasi- vo) region, that is the rectangle with vertices (6 x 10~ 10 ,0.3) 
and (3x 1(T 7 ,3). Indeed, since these models do not predict the overall mass 
scale, we use Am^ to calculate values of the hierarchy factor: 



h 



Am 



assuming that Am^ = 3 x 10~ 3 eV 2 , and then we compare these values of 
h with the calculated ones. The percentage of mass matrices that passes 
the cuts on 623, #13, #12 and h is shown in fig. |l[ As pointed out in ||, |?|], 

2 As usual, M„ = [/*diag[mj x exp(i£j)]U\ rrij < nij+i and Am]j = — mf. In terms 
of the neutrino mixing matrix U (such that vi = UuVi, with i = e,/x,r and i = 1,2,3) 
the mixing angles 9ij, with values 0° < 8ij < 90°, are simply sin #13 = | C C 3 1 , tan 812 = 
\U c2 /U c i\ and tan6» 23 = \U„3/U t3 \. 



2 




10~ 3 10~ 2 10~ 1 10° 10~ 3 10~ 2 10~ 1 10° 
8 8 

Figure 1: Percentage of phenomenologically successful neutrino mass matri- 
ces of the type eq. ([[]) as a function of e. Dashed line denotes the SMA region, 
continuous thin line the LMA region, thick line the low region (last one is 
practically invisible for triplet case). For orientation, we emphasise with 
arrows pointing downward the special cases e = (rrju/m.,-) ' ' ' 2 and with 
those pointing upward the cases e = (sin^c) 1 ' 2 ' 3 ' 4 (#c =Cabibbo angle). 



it is of some importance to distinguish the case denoted as "triplet" in 
the figure, when one generates a matrix with random coefficients in eq. (|l|) 
(namely 0(1) = R\; see [14) for theoretical support) from the case when 



0(1) = R2R3 R%i which corresponds to assume that the light neutrino 
masses are due to the "seesaw" mechanism [15]. In last case, the random 
coefficients are those of the Dirac neutrino couplings, and of the heavy (right 
handed) neutrinos-i?2 and R3 respectively.^] It is seen that a reasonably 
large number of mass matrices of the type eq. (Q) survive the cuts for certain 
values of e. It is rather obvious that for very large values of e the number 
of successful models decreases due to the cut on #13. The lower peak height 
of the SMA- (even more, low-) curves is due to the difficulty to reduce 
sufficiently h (which is alleviated in the "seesaw" case || |7|). In both cases, 
there is a gap with relatively small overlap between the SMA and lma curves. 
This can be explained in the following manner: consider the plane (/t, #12); 

3 In general, both contributions are present: 0(1) = sin 2 77 Ri + cos 2 77 R2Rs 1 R t 2- For 
simplicity, we will discuss only the extreme cases-triplet and seesaw-but it is clear that 
we can bias this or that value of r\ only adding theoretical information. 



3 



Qe 




Qr 


Q e c 




Q T c 


e (degrees) 


3 








3 


2 





.83° 


2 








4 


2 





3.4° 


1 








5 


2 





14.° 



Table 1: Three sets of Froggatt-Nielsen leptonic charges which give an as- 
cending series of e = (v/M)® c values (v/M is discussed in the text). We 
normalise to the lowest charges since we focus only on charged lepton mass 
ratios; Q^c = 2 is needed to reproduce m^/nir; and Q c + Qe c = 6 (or nearby 
values) to reproduce m^/m T = (m e /m r ) 1//3 (correct at the 10 % level). 

as typical values, we have h ~ 1 and 9±2 ~ e, but there is an interesting 
tail (=a less populated region) which shrinks with e where 612 increases 
and h decreases (the correlation being tighter in the "triplet" case). For 
diminishing values of e, the tail first meets the LMA region and only after 
moves toward the SMA region: this creates the gap observed in fig. [I].^] The 
different position and height of the "SMA-peak" from seesaw to triplet case 
indicate just that the distributions (and tails) are different. 

Some models of this type (=values of e) have been discussed already in 
the literature: e = sin 3 "&c hi 01) = "v/ m T in @, 0j = svn n dc hi @! 1 
in |J. Except than in the last case, the stress was put on the correlations 
of the properties of neutrinos with those of charged fermions. However, no 
special emphasis has been put on the case: 

e = $m$c ~ 0.22 or = (m M /m r ) 1/2 « 0.24, (2) 

that we see to be quite interesting in connection with the lma region (in 
both the "seesaw" and the "triplet" case). The case e = m^/rrir w 0.06 was 
emphasised in reference Q; it was shown that the use of U(l) selection rules 
a la Froggatt and Nielsen, with charges Q e = 1, = Q T = for the left 
fields is not in contradiction with the gross features of the charged lepton 
spectrum, if at the same time one assumes Q e c = 2, Q^c = 1 and Q T c = as 
charges for the right fields. The same can be however obtained with other 

4 This is the same argument put forward in sect. 2 of ref. B. But note that eq. (2.1) 
therein is improved assuming that the 33-entry is 2, not 1; at the same time, m v ~ 
(Am^m) 1 ' 2 has to be replaced by m„ ~ (Arriatm) 1 ^ 2 This affects the numerical factor 
in front to eqs. (2.2-4), and the region in figure 1 shifts downward by a factor of 4. 
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choices of charge, as those in table [T] with a different value of the parameter 
that regulates all mixings and hierarchies, v/M = (m^/m T ) ' 5 (note that 
the second model is just the one of Sato and Yanagida j7j, since Q — > 2Q 
but also v/M -> (v/M) 1 / 2 in our table). 

We verified numerically the viability of the alternative choices, and found 
that increasing Q c c by 1 unity also leads to similar results, for slightly larger 
values of e; similar conclusions even decreasing Q e c by 1 unity (but worsening 
a bit the agreement with charged fermion mass ratios in the first two cases of 
table H): Nothing against large values of e as in eq. @ from this side. Indeed, 
it was argued already in |jDif | that the value of v/M should be large, simply 
in connection with the size of the Cabibbo angle; sin $c should presumably 
come out as the result of a fluctuation, if v/M were much smaller. 

Now we tackle one last problem, which is strictly linked to the theoret- 
ical (Froggatt-Nielsen) context. We start with some general consideration: 
One can attempt to classify the models for neutrino masses, depending on 
whether the neutrino mixing matrix 

U = U%U V (3) 

is 1. mostly due to the rotation of the neutral leptons U u , or 2. to the rotation 
of charged leptons Ue, or, finally, 3. if both have a comparable importance 
(at least for some mixing angle). It is easy to argue for the first possibility 
(at least in words); neutrinos look special, charged leptons could always en- 
tail small mixing, and thence U u could play an overwhelming role in U. As an 
example of the second possibility, we quote the "lopsided" models described 
m @. Finally, as an important instance of the last possibility, Ue and U v 
are typically of comparable importance for U in the U(l) approach of Frog- 
gatt and Nielsen fl3]| , since they are controlled by the same charges, namely, 
those of the left leptonic fields (note, however, that a partial degeneracy in 
the neutrino mass matrix might enhance the role of U v ). Strictly speaking, 
most of the considerations above apply to the first case, in that (following 
previous studies) we simply ignored the role of the mixing due to the charged 
leptons (Ue ~ 1). However, since the approach with random numbers is 
naturally (though not unavoidably) connected with the use of U(l) selection 
rules for fermion mass matrices, and since the existing studies of neutrino 
mass matrices of this type @, ^, H |l0| do not take into account the point, we 
decided to investigate what is the effect of Ue in models with U(l) selection 
rules. Thence, we additionally generated the random mass matrices Rq as 
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Table 2: Calculated neutrino properties for the U(l) models of tab. |l[ in 
the cases with triplet or seesaw (t and s resp.) mass mechanism, and with 
or without the account of the lepton mixing matrix Ue (w and w/o resp.). 
The 3 parts of the table correspond to the models defined in tab. |l| (in the 
left-upper corners, the values of £ in degrees). All angles in the table are in 
degrees, while SMA, LMA and LOW denote percentages (cuts as in fig. |l]). 



coefficients of the charged lepton mass matrix = U E diag[m^] Ve '■ 



W/Ie oc diag 



M 



0(1) diag 



AI 



(where 



e, (jl, t) enforcing the values of the U(l) charges in tab. [I] to 



specify the models fully. Only if the mass ratios (m e /m T ) and (m^/m T ) are 
reproduced within 30 %, we calculate Ue and estimate its effect on U in eq. 
(||) (it is important to implement such a condition to gauge out cases when 
the mixing angles come out artificially large). 

The results are presented in table ^. The spread^] in 623 does not change 
much with e, however, it increases significantly with the inclusion of Ue 



5 The spread 6x on the quantity x is calculated as 8x 2 = {x% — (x)) 2 /N, where Xi 
is the result of the i th simulation, (a;) = Xi/N is the average, and iV = 10 6 . 
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effects; it is a pity that these simple models are unable to give an indication 
on the size of the deviation of #23 from 45°, one of the most interesting 
quantities to be searched in future experiments. The average values of 9%2 
and also increase with the inclusion of Ue- If Ue effects are not included, 
#13 is on average smaller than e [|7j ; the reason is just the numerical factor ~ 2 
discussed in footnote ||. The angle #13 is rather close to e when these effects 
are included, and close to the experimental limit when e is large. Note that 
even a seemingly modest increase in #13 , say by a factor of 2 (table ||, triplet 
case) is an important message for the searches of v e appearance in terrestrial 
experiments, since the probabilities of conversion in vacuum depend strongly 
on this parameter: Piy^ — > u e ) cx (#i3) 2 . In tab. ^ we also present the value 
of the following quantity: 

|(M„)ee| 

mee " (Am§! ) 1/2 

These calculations show that, for the larger value of e considered, the ee- 
element of the mass matrix |(M 1/ ) ee | can reach the several meV level. This 
is tantalisingly close to the expected sensitivity of the next generation neu- 
trinoless double beta decay experiments, namely 10 — 20 meV. 

Let us summarise and discuss the results. Neutrino masses are expected 
to be non-zero, and indeed the atmospheric neutrino mass scale ~ 50 — 60 
meV is very well compatible with the ideas of grand unification, but the size 
of the mixing angles is a puzzle, in particular the strong indication that at 
least one of them is large-maybe maximal. We discussed a class of neutrino 
mass models (eq. (|])) that are inspired to the principle that all elements of 
the mass matrices (seen as fundamental quantities) should be equally large, 
unless explicitly suppressed. We emphasised the type of models that are 
phenomenologically successful, and discussed the stability in the choice of 
coefficients by the help of random number generators (fig. [[]). We discussed 
also the connection with known properties of the charged fermions (and, in 
particular, leptons), and remarked that large values of the "order parameter" 
e (see eq. (§)) are quite interesting in connection with the LMA solution 
of the solar neutrino problem. In present context, the triplet mechanism 
is not disfavored in comparison with the seesaw mechanism for neutrino 
mass generation; one could argue instead that the triplet mechanism is more 
predictive, for the simple reason that it favors solar neutrino solutions with 
as large h as possible. In close observance to the Proggatt-Nielsen approach, 
we also pointed out the relevance of the mixing due to charged leptons for 
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the neutrino mixing; this effect unfortunately renders almost random the 
expected value of the parameter #23 in these models, though the maximal 
mixing #23 = 45° is the most likely possibility. Large values of the mixing 
angle #13, and perhaps an observable |(M I/ ) ee | can be naturally incorporated 
in the model, if the parameter e is large as in eq. (g). 



I thank W. Buchmiiller, D. Delepine, A. Grillo, P. Lipari, H. Nunokawa, 
G. Senjanovic, M. Tanimoto and especially Z.G. Berezhiani for clarifying 
discussions. After these calculations were completed, we learned of an in- 
teresting study [17] of illustrative cases of the mixing matrices Ue and U u ; 
when comparison is pertinent, our conclusions agree with theirs. 



Appendix: What is a Random Coefficient? 



In the spirit of the approach, the random coefficients of the mass matri- 
ces should not be enumerated among the parameters; however, their form 
has to be fixed in order to perform the scan of possible models. To give 
an idea of the effect of changing the generator, we compare the percentage 
of successful models in (sMA,LMA,LOW)-regions for the model of Sato and 
Yanagida (seesaw case, Ue included-tab. [| middle part, last line). Varia- 
tions amounts to a factor of a few in most extreme cases: (0.24, 0.16, 0.11) % 
for Z , (0.42,0.22,0.11) % for Z x . We denote with Z& the numbers in the 
complex plane with random phase, and with modulus uniformly distributed 
in an interval [—5, 5] around 1 [Z\ is a circle in the complex plane, Zq the 
circumference). Following Q, we adopt in our study the choice Z$, with 
5 = 0.2, clearly consistent with the view that the coefficients are "order 
unity" . 

There is another point of ambiguity: How to treat a symmetric random 
mass matrix R (indeed, the matrices R\ and i?3 described in the main text 
are symmetric). One possibility is to generate just the elements Rij with 
i > j, and then set Rji = Rij] another one is to generate the full matrix, and 
then replace Rij and Rji by their average.^] Following the previous studies, 
we adopted the first prescription, and noted that the second prescription 

6 These are clearly different prescriptions; for instance, the sum of 2 numbers uniformly 
distributed is not uniformly distributed. 
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leads to results that differs little from those in tab. Q (in the same case of 
previous paragraph); 623 is the same; #13 = (3.2 ± 1.6)°, 6*12 = (5.6 ± 5.0)°, 
h = .11± .15, m ee = (4.7±3.9) x 1CT 3 ; finally, the sma,lma and low success 
percentages are 0.26,0.10 and 0.08 %. 

Now, a delicate aspect; in this work (following previous studies) we as- 
sumed that all coefficients 0(1) are distributed in the same manner, but a 
priori it is unclear whether this assumption is fair. Indeed, if terms higher 
order in v/M arise through the exchange of virtual particles as suggested in 



[13|, there might be not only a piling up of (v/M)'s but also of the coeffi- 
cients themselves]^] We do not want to deny the interest of this simplifying 
assumption (that we adopt in this study), but we point out a risk of minor 
reliability of the predictions which depend essentially on higher order terms. 
From this point of view, the prediction of a relatively large value of h and 
of 623 are more reliable than those on the other mixing angles; the one on 
|(M^) ee | has the greatest dependence on this crucial assumption, and might 
be considered, then, even less reliable. 

A final warning, of more general nature: Mass matrices with random 
coefficients can help to emphasise certain possibilities that are compatible 
with present (lack of) information, but one should be careful to interpret 
the results of this type of calculations in "probabilistic" term. Indeed, if, 
for consistence, the absolute scale of the neutrino mass matrix was also let 
fluctuate; or also, if one required that the (masses of the) charged leptons 
were reproduced within experimental errors; etc.; the probability of success 
of these attempts would have been practically zero. Similarly, if a theory of 
the coefficients order unity were given, any "statistical" consideration (like 
the present ones) would have been much less relevant. 
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